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Abstract 

We analyze two kinds of matched filters for data output of a spherical resonant GW detector. 
In order to filter the data of a real sphere, a strategy is proposed, firstly using an omnidirectional 
in-line filter, which is supposed to select periodograms with excitations, secondly by performing 
a directional filter on such selected periodograms, finding the wave arrival time, direction and 
polarization. We point out that, as the analytical simplifications occurring in the ideal 6 transducers 
TIGA sphere do not hold for a real sphere, using a 5 transducers configuration could be a more 
convenient choice. 
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I. INTRODUCTION 



Spherical resonant gravitational wave (GW) detectors are the evolution of bar resonators. 
For an equivalent size they have a higher-cross section than other geometries and are the 
only kind of GW detectors able to to determine the signal arrival direction and polarization, 
even operated at a single location on Earth. Nowadays two spheres are under development 
in the Netherlands (MiniGRAIL)fl and in Brasil (MARIO SCHENBERG)^. Their exper- 
imental situation at present requires an efficient and concrete data-analysis strategy, which 
supplement the available literature (see for exemple Works in this direction have 

already been done, see for example [5(. 

In this paper we will focus on the case of matched filtering for such a detector, starting 
from some basic notions of Id]. Matched filtering has been shown to be a key point in 

q n 

the burst oriented data analysis for resonant detectors [7|] and for interferometers [8[. We 
will discuss the question from a theoretical point of view, keeping in mind the technical 
implementation of our propositions. We propose a strategy in order to filter the sphere data, 
firstly by an omnidirectional in-line filter, which we call "refined energy matched-filter" , in 
order to select signal excitations, then by performing a "directional matched-filter" on events 
selected by the first filter, extracting the wave arrival time, direction and polarization. 

In actual experiments, the preferred configuration for the transducer placement is the 
truncated ichosaedral (TI) proposed by Merkowitz and Johnson jjj. We find that when 
data analysis is performed for a realistic sphere with resonators and electronic readout, the 
analytical simplifications of an ideal sphere with transducers placed in the TI configuration 
do not hold anymore. Therefore using a 5 transducer configuration as the pentagonal hexa- 
contahedron (PHC) proposed by Lobo and Serrano [h]] can be a more convenient possibility, 
because simpler to build and manage. 

This paper is structured as follows: in the next section we introduce notations and a 
schematic model of spherical detector. In section III II and [IV] we present two different ways 
of implementing a matched filter with a description of advantages and limitations for each 
one. In the last section we draw a strategy in order to perform data analysis with these two 
filters and end with our main conclusions. 
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II. MODEL OF SPHERICAL GW DETECTOR 



Analytical considerations 



We will provide a general frame to model a GW spherical detector equipped with me- 
chanical resonators and electronic readouts. We only take into account couplings which are 
linear in each of the considered variables. 

In accordance with previous literature, we call z the amplitudes of the five fundamental 
normal modes (z is then a five- component vector) and q the positions of the six mechanical 



resonators in the TI configuration 9|, . The model is completed by the current readouts 
I, which is a vector of k ■ 6 components, with k = 2 for a system of capacitive single-SQUID 
transducers 



In Fourier space the equations of motion for all those variables take the matrix form 

HQ 



(i) 



where F$ and represent the noise forces acting on the sphere and the resonators, V the 
noise voltages and Fqw the GW effect. We use this equation in the form 
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We have work with transducers in the TI configuration, the most used in recent literature, 
but this is not a unique choice. In this special case one can indeed easily find out the so 
called "mode-channels" through the formula 



p = BI 



out; 



(3) 



lout being the 6-vector containing the last six components of I, and B the 5x6 matrix 
representing the transducer positions. Its explicit expression can be found, for example, in 



14i | . In the following we will always work with the output currents and therefore drop the 



subscript "out". 



3 



B. Numerical model 



The numerical model used in this paper is a direct sequel of those developed in 12j. In the 
following we will restrict ourself to a 1 meter radius CuAl filled sphere, with a quality factor 
for the sphere's fundamental modes Q$ = 5 • 10 7 with 6 capacitive transducers having both 
mechanical and electrical modes tuned on the sphere mode. We have tuned the transducers 
to the quantum limit (A^ phonon = 1) rather than to use the present best sensitivity value 
(-Wphonon ~ 50) but this does not affect our analysis. 



III. DIRECTIONAL MATCHED FILTER (DMF) 



In a noise-free detector the output response to a GW burst is given in Fourier space 
by Q 



I 



-l/2uo 2 m s R sX GT v {9,(t ) ^) 



(4) 



where m s , R s \ are the mass and the effective length of the sphere, G the 6x5 sub-matrix 

fh+\ 

of G relating the forces onto the modes to the output currents and Ty(#,0, ip) _ the 

\h x J 

projection of the GW tensor on the different modes, which depends on the wave polarization, 
ip, and arrival direction, (6,<f)). Choosing a coordinate system with origin on the center of 
the detector and the z axis pointing the zenith, 9 is defined to be the angle between the z 
axis and the wave direction and <fr the angle between its projection on the xy plane and the 
x axis. Using the tensor spherical harmonics formalism, which allows a definition for the 



five fundamental modes 



13], Ty can be written as 
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We can then perform the matched filter searching the expected pattern, say Igi, into I, 
the output of the detector. This is performed using the noise weighted inner product 6] 



(ii|i 2 ) = / ^iiHS-V)^ 



(6) 
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where S, the noise matrix, is the expectation value of the product I no iseI no ise with I no i S e 
representing the output due to noises only [6]. For experimental purposes, the noise matrix 
can be computed as an average over the detector output S = (11^). The output due to 
noises is given by 



fF s \ 



I, 



G 



(7) 



fa 

V v/ 

with G being the sub-matrix given by the last 6 lines of G. In order to rewrite the matrix 
S in terms of G and the noise spectral densities, we use the statistical properties and the 
independence of the noises F5 f# and V [12 ] . 
We obtain the filtered signal as 

I^aOS-WCw) 



7(0;) 
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(8) 



Ifli is the filter output obtained by choosing a waveform and a polarization of h+ l (t) , h^\t) 
and a direction for the expected wave. The coefficient of the projection on the pattern is 

(lm|l) 



, dw , s 
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The same equations hold using mode- channels: one has just to multiply eq.(j4|) on the left 
by B and perform the replacement I — ► p. 



A. Limitations 



The principal problem with this kind of matched filter is that pm is a function of the 
arrival direction and of the polarization angle. Therefore such a filter can be built for a 
specific GW direction, polarization and shape (usually matched filters are built for 5-like 
excitations and we will perform our test with such a waveform but one can also use sin- 
Gaussian or a specific template) which therefore are supposed to be known "a priori" . This 
is not a realistic possibility. Such a filter is unable to detect with high efficiency an excitation 
with a different arrival direction or polarization. 

The loss of sensitivity of the directional filter constructed for the direction 9, <p and 
polarization ip when a GW is coming from another direction, say 9' and 0', with another 
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FIG. 1: Direction dependence of s on the cylindrical coordinates 9' , (/>' . 9' is the radial coordinate, and <fi' 
the angular coordinate. We have fixed 9 = 0, = 0,-0 = 7r/4, ft.?} 1 = 1, = and plot s for 6*' e [0,7r/2] 
and 4>' € [0, 2ir]. In order to avoid a mulitvalued function at 0' = we have corrected the polarization angle 



12f. 



setting tp' = ip + (j)' . We use the sphere model with capacitive transducers described in 

polarization, ip', is obtained by fixing h and computing the ratio 

s — _ — (10) 

(i ffl (M,V0|ifii(M,V0) 

In Figure [U, we plot the 9' and 0' dependence of s for our sphere model. 

The polarization dependence of s, for a given direction (9 = 9', <fi = </>'), can be computed 
using the explicit ip dependence of eq.((5]). Without loss of generality we can fix = 1, 

= and ip' = ip + Sip, getting 

/ cos 2ib sin2^\ t / cos(2ib') sm(2ib') \ /1\ 
s = (1 0) M 

V-sin2^ cos 2^/ \ - sin(2^') cos(2^') / \0 J 

= - (sin(4^ + 25ip) + sin(2^)) M n + 

- (- cos(4^ + 25^) + cos(2<ty>)) M 22 . (11) 

All terms which are independent of ip are put into the 2x2 matrix M. As an exemple, for 
ip = 7r/4 the sensitivity of the filter drops as cos(25ip). For each direction, reconstructing 



the polarization is possible by fixing a polarization angle ip and performing two identical 
filters up to a shift of 7r/4 in the polarization angle. 

If the wave arrival direction is not known a priori (as it is always the case), one possibility 
is to produce several filters for different directions. The detector output has then to be 
processed with all these filters, giving finally a direction and a polarization which maximize 
r. The clear advantage is that this procedure gives an estimate of both the arrival direction 
(the one corresponding to the "better" filter) and the polarization. However this procedure 
needs high computing power since it requires 6 x N more computational time than resonant 
bar matched filtering, N being the number of filters. We can estimate N using Figure [U 
If we want an efficiency of say s > 0.8, we can deduce that a filter with arrival direction 
in 9 = can only see signals from directions with 6 < n/6 and therefore cover ~ 1/8 of 
the hemisphere (remember that antipodal directions are identified). Taking into account the 
polarization, we obtain that N has to be at least of order 20. However, once a burst has been 
caught by the filter, a hierarchical algorithm has to set up (trying maybe also some other 
wave forms) in order to reconstruct an accurate arrival direction. This is not satisfactory for 
a in-line analysis but can be useful in order to analyze particularly some interesting events. 



B. Performance 



The sensitivity of a detector depends on the method used to analyze its data. We quantify 
now the sensitivity that can be reached with the directional matched filter by computing its 
"strain sensitivity", defined as the amplitude h c (u) that an incoming GW needs in order to 
produce the same output as detector internal noise, that is to have a signal to noise ratio 
(SNR ) of 1. As expected, for matched filtering the highest SNR is reached when the filter 
has the same form as the signal (when I is equal to Igi), see js] and references therein. This 
situation corresponds to optimal filtering. In this case we can define the filter SNR density 
as 

a(u) = ^(^S^HI+H . (12) 

Furthermore, without loss of generality, we can perform a rotation of the polarization frame 
in order to set one of the polarizations equal to 0, 

i> - V' (13) 
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FIG. 2: Strain sensitivity computed for directional filter. We use the l[m] radius sphere model with 



capacitive transducers of 



12 ] and display the different noise contributions. We have used a r5-hTter 







(14) 
(15) 



Using these new variables together with eq. (j4]), we get 

= f{ U )\h a {uj)\ 2 (16) 

where f(u) is a proportionality coefficient. Setting (j{uj) = 1 gives the strain sensitivity 

/ 1 

h c {u) - 



(17) 



In Figure El h c (u) is plotted for a directional matched filter applied to a l[m] radius sphere 
with six single SQUID capacitive transducers in TI configuration as described in 12]. We 
stress that such a strain sensitivity is reached only for optimal filtering, and it will be reduced 
if the filter is not built with exactly the same arrival direction and waveform as the received 
GW signal. 



IV. ENERGY MATCHED FILTERING 

We observe that the sum of the squared excitations of the five fundamental sphere modes 



(h* + h* x ) T V T 



h., 



(18) 



is proportional to \h + \ 2 + \h x \ 2 and independent of the direction. Furthermore, one of the 
most useful properties of the mode-channels formalism is that the vector p is proportional 
to the fundamental mode excitations of the sphere, at least when considering only an ideal 
sphere with TI resonators. In other words, this means that BG is proportional to the 



identity. The proportionality factor can indeed be easily worked out [141 ]. 

As a consequence G^B^BG is proportional to the identity too, and the norm of the vector 
p in Fourier space 



fhf 

■v 
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(19) 



has the same kind of proportionality to \h^\ 2 + \h^\ 2 and directional independence 151 ]. 

The main idea of this section is to use this directional independence in order to perform 
a matched filtering on the total energy stored into the sphere, which can be computed from 
the detector signal 

E(u) = p f p = I+B+BI . (20) 
The noise weighted inner product is defined as 

El \E 2 ) = [ ^(fE^) 

where the noise contribution 

-Enoisc = Inoise-E^-^Inoise (*^2) 

can also be computed by means of the noise spectral densities. 
For a given signal E^, the filter gain is then 



and the filter is given by 



Integrating we get 



R=(E m \E m ) (23) 



1E[U) = R^IJJV- (24) 



noise \ 
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We define the efficiency with respect to the directions and polarizations in the same way 
as we have done for the directional filter: 

_ (£ m (M,V0|£m(fl',4>',^)) 

(^(M,V0|£*i(M,V0) ' 

Since we have defined the energy matched filter to avoid angular dependence, we expect 
that se = 1. 

Once the filter waveform has been fixed, its amplitude is determined by a single parameter 
(h m ) 2 . Typically, for a 5- like excitation we have |#l(o>)| 2 + |^(o;)| 2 = (h m ) 2 = constant and 
for a sin-Gaussian \hf(t)\ 2 + \h^(t)\ 2 = (h m ) 2 sin 2 (u;ot)e- 27(t - to)2 . We can build the filter 
Efu(u) from eq. ( TT9i) with the chosen pattern and h m = 1. In this case, computing Te on 
a detector output periodogram E returns the amplitude of an eventual hosted excitation 
having the same waveform as the filter: 

T E = h 2 . (27) 



A. Limitations 



The energy matched filter is efficient to determine the presence of an excitation out of 
noise, independently on the direction and polarization of the wave. Consequently it cannot 
be used in order to determine the arrival direction of the GW but rather as a in-line filter 
able to detect efficiently the presence of a GW signal. However, this supposes that BG is 
proportional to the identity, which is not true for a real spherical detector. For a physical 
model of resonant spheres, we have to include the splitting in frequency due to the Earth 
gravity field and electronic readout. 

It is interesting to split BG in order to quantify the deviation to a multiple of the identity 

N(w) = BG(w) - C(u)l 5 , C(u) = ^ y £(BG(u)). i (28) 

i=l 

where I5 is the 5x5 identity matrix. 

Numerical investigations for models with capacity transducers, show that the elements 
of N may be non negligible against C. In particular for our model the magnitude of the 
biggest elements of the matrix N is typically of the order of 7% of C, but can reach 400% of 
C for uo near the resonances. As a sequel, E(u) is direction dependent, as shown in Figure 

El 
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FIG. 3: Under the same conditions than Figure[T]we plot the direction dependence of ^/se- Note that y/sE 
vary only between ~ 1/2 and 1. 



B. Refined energy matched filter (REMF) 

We adapt our analysis in order to recover the direction-independence of E(u). With this 
aim, we change the mode-channel definition ([3]). We need to replace the matrix B by a 
matrix B(u;) such that the product B(o;)G(o;) is proportional to the identity. We define the 
matrix / 1 \ 

1 

B b 1 5 



V 



V 



and the new mode-channels p(k>) are defined as 



d a 



(29) 



J 



The choice of the last column is arbitrary. For each value of uo we can compute the 5x5 
matrix b as ( 1\ 

1 
1 
1 



GM 



/ 



(30) 



p(u) = B(w)I(w). 



(31) 
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Note that the matched filter does work since we fix a unique set of matrices B(u;) replacing 
the matrix B in equations from (fT9l) to ff27|) . 

The main consequences of this change are the following: 

• We need the knowledge of G(a>) in order to compute the mode-channels. 

• The mode-channels are only defined through the Fourier space computation. 

• The product B(u;)G(u;) is proportional to the identity but the proportionality constant 
is now C(u), an arbitrary (but non zero) function of uj. 

• Je(w) has to be manipulated with some caution since it depends on C(u). However 
the property Ye = h 2 remains true and is sufficient in order to select periodograms 
hosting excitations. 



C. Performance 



We checked with our numerical model and for 5-like waveform that the redefined matched 
energy filter is able to give a response which is independent of the polarization and the 
direction (of both the filter function and the incoming wave) with relative error compatible 
with numerical errors O(10~ 16 ). 

We expect the strain sensitivity corresponding to REMF to be worse than the one ob- 
tained from DMF. In particular, the computation of E as a sum of the squared mode- 
channels, implies that we loose the relative phase and therefore the possibility of performing 

m 

a coherent analysis [12j. In order to compute the strain sensitivity for the REMF, we need 
the SNR density <7e(u;), which is given by the integrated term in eq. (1231) . This is, by con- 
struction, proportional to (j/i+| 2 + l^x'l 2 ) an< ^ the strain sensitivity is then obtained as the 
4th- root of the proportionality coefficient. Figure H] shows the strain sensitivity of the refined 
energy matched filter. Note the presence of horns, typical of non-coherent analysis 12] . 

We emphasize that even if the sensitivity is worse for the energy filter, this is partially 
compensated by the fact that the energy filter performance is the same for all the directions 
and polarizations. 
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FIG. 4: Strain sensitivity obtained for the energy filter. This computation is performed for the same sphere 
model as in Fig. [3] 

V. CONCLUSIONS AND FUTURE DEVELOPMENTS 

In this paper we have analyzed two ways for matched-filtering data output of a spherical 
GW detector. In section [IV] the "refined energy matched filter" has been presented, whose 
efficiency does not depend on the GW direction. 

For the specific setup of each existing GW spherical detector, it would be worth to inves- 
tigate in more details the performance of the pipeline data analysis proposed in this work: 
firstly using the REMF as an in-line filter in order to select periodograms with excitations. 
REMF requires some computer memory, since the set of matrices B(u;) has to be stored, 
but only few more stages of computation than the matched filter used for the cylindrical 
detectors. 

Then, on periodograms selected by REMF, the directional matched filter presented in 
section III 1 1 of this paper can be used in order to extract the arrival time, direction and 
polarization of the GW signal. Furthermore, once we have obtained (for a chosen ip) h + 
and h x from the DMF procedure, we have a check of the coherence of the whole procedure 
by making the comparison with the value of \h + \ 2 + \h x \ 2 obtained from the refined energy 
matched filter. 

Finally we stress that the mode-channel simplifications occurring when the 6-transducers 
TI configuration is considered do not hold in a realistic model of a spherical antenna. 
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Therefore, for the data analysis strategy proposed in this work, it may be convenient to 
choose working with only five transducers. In such a case the G is a 5 x 5 matrix and, 
finding a configuration of the transducer such that it is invertible for all values of u, we can 
perform a REMF with B = Gr -1 . 

We would like to thank C. A. Costa and A. Malaspinas. This work is partially supported 
by the Swiss National Science Foundation (FNS). 
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